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This  article  is  a  literature  review  of  use  of  computational  fluid  dynamics  (CFD)  codes  to  model  the  in¬ 
cylinder  fluid  flow,  turbulence  and  spray  characteristics.  This  study  is  based  on  the  reports  of  about  60 
scientists,  who  published  their  results  between  1978  and  2008.  Most  of  the  scientists  and  researchers 
used  CFD  codes  to  analyze  the  models  under  simulation  conditions  and  compared  these  simulated  results 
with  experimental  results.  Some  scientists  reported  that  different  engines  exhibit  different  behaviors 
with  similar  fuel  sprays  and  Re-Normalized  Group  (RNG)  k-s  model  is  the  best  applicable  turbulence 
model  for  engine  simulation.  The  KIVA  code  is  widely  used  for  model  development  in  academia  due  to  the 
availability  of  the  source.  However,  its  capability  for  resolving  complex  geometries  is  limited.  On  the 
other  hand,  other  commercial  CFD  codes  such  as  STAR-CD,  FIRE,  VECTIS  and  FLUENT  are  frequently  used 
by  the  industry  due  to  their  superior  mesh  generation  interfaces  and  because  of  their  available  user 
support.  Some  scientists  combined  STAR-CD  and  KIVA  code  for  the  engine  simulations  but  they 
concluded  that,  it  would  be  preferable  to  implement  the  advanced  submodels  directly  into  one 
commercial  code  for  engine  simulations. 
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1.  Introduction 

Fluid  flow  in  an  internal  combustion  engine  presents  one  of  the 
most  challenging  fluid  dynamics  problems  to  model.  It  is  because  the 
flow  is  associated  with  large  density  variations.  Further,  the  fluid 
motion  inside  a  cylinder  is  turbulent,  unsteady,  cyclic  and  non¬ 
stationary  both  spatially  and  temporally.  The  combustion  char¬ 
acteristics  were  greatly  influenced  by  the  details  of  the  fuel 
preparation  and  the  distribution  of  fuel  in  the  engine  cylinder, 
which  was  mainly  controlled  by  the  in-cylinder  fluid  dynamics. 
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Further,  fuel  injection  introduces  additional  complexities.  It  was 
mainly  due  to  the  difficulties  in  understanding  the  physics  of  two- 
phase  flows  that  vary  spatially  and  temporally  from  very  dense  to 
fairly  dilute.  Pollutant  emissions  were  controlled  by  the  details  of  the 
turbulent  fuel-air  mixing  and  combustion  processes.  So,  a  detailed 
understanding  of  these  processes  is  required  to  improve  perfor¬ 
mance  and  reduce  emissions  without  compromising  fuel  economy. 

2.  In-cylinder  fluid  flow  models 

The  role  of  in-cylinder  air  motion  commences  from  the  very 
start  of  the  engine  cycle.  During  the  intake  stroke,  the  incoming  air 
generates  flow  structures  with  large-scale  turbulent  motions 
within  the  cylinder,  which  in  turn  determines  the  extent  of  mixing 


S.A.  Basha,  K.  Raja  Copal / Renewable  and  Sustainable  Energy  Reviews  13  (2009)  1620-1627 


1621 


between  the  fresh  charge  and  the  residuals,  as  well  as  internal  and 
external  heat  transfer  rates.  The  flow  field  plays  a  vital  role  in 
charge  preparation  and  conditioning  for  ignition.  During  the 
combustion  processes,  the  flow  field  continues  to  exert  its  effect  on 
heat  transfer,  peak  flame  temperatures  and  thermal  stress  levels.  In 
the  post  combustion  period,  the  nature  of  fluid  motion  determines 
the  level  of  pollutant  emissions.  The  flow  field  also  determines  the 
cyclic  variation  exhibited  by  it  and  all  other  processes,  which  it 
influences. 

Margery  investigated  the  affect  of  intake  duct  length  on 
volumetric  efficiency  and  the  in-cylinder  flow  field  in  a  four- 
stroke  single  cylinder  direct  injection  (Dl)  diesel  engine  motored  at 
1000-3000  rpm.  Three  lengths  of  straight  duct  upstream  of  the 
helical  inlet  port  of  the  engine  are  considered  and  measurements  of 
instantaneous  mass  flow  rate  and  pressure  drop  across  the  port  are 
reported  as  a  function  of  duct  length  and  engine  speed  [1],  Zhang 
et  al.  analyzed  the  combustion  phenomenon  of  a  direct  injection 
diesel  engine  by  using  the  cross-correlation  method  and  two-color 
method  to  measure  the  combustion  flame  motion  and  the  flame 
temperature,  respectively  by  processing  the  high-speed  in¬ 
cylinder  photographs.  The  effects  of  engine  parameters  such  as 
pumping  rate,  injector  nozzle  hole  size  and  injection  timing  on 
combustion  processes,  particularly  on  flame  motion  and  flame 
temperature  were  studied  [2]. 

Chen  and  Dent  presented  a  combined  experimental  and 
computational  study  of  steady  flow  through  a  curved  inlet  port. 
Flow  structures  affected  by  valve  lift  and  port  shape  are  predicted 
by  the  computational  fluid  dynamics  (CFD)  code  STAR-CD.  Local 
static  surface  pressure  maps  are  obtained  experimentally  [3].  Kong 
et  al.  developed  an  image  acquisition  and  processing  camera 
system  to  acquire  luminous  flame  images  from  the  combustion 
chamber  of  a  single  cylinder  heavy-duty  diesel  engine.  Experi¬ 
ments  were  conducted  both  single  and  split  injection  cases  at  90 
and  45  MPa  injection  pressures  at  3/4  load  and  1600  rpm  with 
simulated  turbo  charging.  The  experimental  results  are  compared 
to  numerical  computations  using  a  version  of  KIVA-II  code  with 
improved  spray,  ignition,  combustion  and  emission  models.  Good 
levels  of  agreement  are  obtained  with  measured  cylinder 
pressures,  heat  release  rate,  ignition  timings  and  locations,  and 
flame  shapes  [4j. 

Zhang  et  al.  carried  out  an  experimental  study  aiming  to 
investigate  the  effects  of  combustion  chamber  geometry  on 
combustion  process  in  an  optically  accessible  direct  injection 
diesel  engine.  Flame  movement  behaviors  such  as  the  distribution 
of  flame  velocity  vectors  and  the  averaged  flame  velocity  inside 
and  outside  the  combustion  chamber  are  measured  by  means  of 
cross-correlation  method.  The  combustion  geometries  considered 
were  dish  and  re-entrant  chambers  [5].  Stephenson  and  Rutland 
assessed  the  differences  in  the  in-cylinder  flow  field  present  during 
combustion  resulting  from  different  intake  flow  configurations  and 
compared  with  the  significance  of  spray-wall  interaction  effects 
on  combustion  and  emission  performance.  The  key  differences 
observed  in  the  flow  field  at  the  onset  of  fuel  injection  were 
turbulent  kinetic  energy  and  the  magnitude  of  the  swirl  flow. 
Quantities  such  as  the  tumble  and  swirl  flows  and  crossing 
frequency  length  scales  proved  to  be  helpful  in  understanding  the 
development  of  the  flow  field  during  intake  and  compression.  By 
the  time,  the  fuel  injection  began,  most  of  the  tumble  flow  angular 
momentum  had  been  dissipated  or  converted  into  swirl  flow.  The 
most  significant  differences  in  the  injection  flow  field  were  in  the 
turbulence  intensity  and  swirling  flow,  both  of  which  maintained 
differences  in  intensity  developed  early  in  the  intake  flow 
throughout  the  intake  and  compression  strokes.  At  idling 
operation,  the  differences  generated  by  intake  flow  were  found 
to  be  considerably  less  important  than  the  effect  of  fuel  spray-wall 


impingement  in  respect  of  pollutant  formation.  At  75%  load,  the 
effects  of  intake  flow  were  comparable  to  those  of  wall  interaction 
[6], 

Stephenson  et  al.  conducted  a  parametric  study  on  the  effects  of 
swirl,  initial  turbulence,  oxygen  concentration  and  ignition  delay 
on  fuel  vaporization,  mixing  and  combustion  process.  The  effects  of 
intake  flow  field  was  further  investigated  by  varying  the  geometry 
of  the  intake  ports  and  intake  runners  in  a  model  dual  port  direct 
injection  diesel  engine  in  KIVA  simulations.  The  code  solves  the 
three  dimension  conservation  equations  and  uses  an  improved 
RNG  k-e  turbulence  model  [7],  Catania  and  Spessa  used  an 
advanced  hot  wire  anemometer  (HWA)  technique  for  air  velocity 
measurements  in  the  combustion  chamber  of  an  automotive  diesel 
engine  with  a  slightly  re-entrant  cone  shaped  in-piston  bowl  and 
one  helicoidal  intake  duct  under  motored  conditions.  The  in¬ 
cylinder  flow  data  are  acquired  in  a  crank  angle  interval  ranging 
from  30  degrees  before  intake  top  dead  centre  (TDC)  to  60  degrees 
after  compression  TDC,  so  as  to  investigate  the  direct  and  reverse 
squish  flows  produced  by  the  bowl,  an  intense  and  high  frequency 
turbulence  production  occurred  during  the  induction  stroke  which 
was  ascribed  to  the  secondary  and  main  intake  jet  flows.  Over  the 
end  of  this  stroke  and  during  the  start  of  compression,  the 
turbulence  intensity  and  frequency,  the  mean  flow  and  its  cyclic 
fluctuation  showed  to  decay  at  crank  angles  and  rates  that 
increased  with  the  engine  speed.  During  the  main  part  of 
compression,  these  quantities  tended  towards  almost  stationery 
distributions  with  a  fluctuating  pattern  at  the  highest  speeds  and 
with  a  maximum  in  the  mean  flow  at  the  lowest  speeds  over  mid 
compression.  These  trends  were  most  likely  due  to  flow 
unsteadiness,  compression  and  stratification  effects.  At  the  end 
of  compression,  a  further  decay  occurred  in  the  mean  flow  and 
turbulence  intensities  up  to  the  onset  of  the  direct  squish  flow 
where  an  increase  of  intensity  and  frequency  are  generally 
observed.  During  the  investigated  part  of  the  expansion  stroke, 
a  strong  growth  in  the  mean  flow  and  in  both  turbulence  intensity 
and  frequency  are  detected  at  all  speeds.  This  is  ascribed  to  the 
conical  in-piston  bowl  nozzle  effect  of  intensifying  the  reverse 
squish  jet  flow  and  its  interaction  with  swirl  [8]. 

Leylek  et  al.  examined  steady-stake  intake  region  flow  for  low, 
medium  and  high  valve  lift  cases  in  a  caterpillar  diesel  engine,  with 
the  intent  of  investigating  total  pressure  losses.  Computationally 
predicted  loss  pockets  are  identified  and  presented  thoroughly  for 
high  valve  lift  case.  Result  shows  upto  30%  of  total  pressure  loss 
occurring  in  the  upstream  of  valve  clearance  region  and  highest 
losses  occurring  in  the  valve  clearance  region  [9], 

Bo  et  al.  assessed  the  accuracy  of  prediction  of  the  air  motion  and 
spray  characteristics  in  a  direct  injection  diesel  engine  with  the  CFD 
code.  The  Laser  Doppler  Anemometiy  (LDA)  measurements  of 
velocity  field  are  obtained  for  the  induction  and  compression  strokes 
upto  the  time  of  ignition.  The  predictions  of  gas  velocities  spray 
penetration  and  droplet  velocities  for  the  most  part  reproduce  the 
experimental  data.  However,  further  evaluation  of  the  spray 
modeling  particularly  in  respect  of  the  effects  on  fuel  vapor  and 
gas  temperature  needed  [10].  Chen  et  al.  simulated  the  transient 
flow  for  the  induction  stroke  using  STAR-CD  version  2.21  software. 
The  standard  k-e  model  was  employed  to  represent  the  fully 
turbulent  flow  with  the  ‘wall  function'  for  description  of  boundary 
layer  behavior  [11].  Auriemma  et  al.  made  in-cylinder  measure¬ 
ments  under  motoring  conditions  in  a  light  duty  diesel  engine 
equipped  with  a  re-entrant  bowl  in  piston  combustion  chamber. 
Tangential  and  radial  components  of  the  air  velocity  are  acquired 
over  a  crank  angle  range  of  90“  above  TDC.  The  mean  motion,  integral 
time  scale  and  Reynolds  shear  stress  are  estimated  by  applying  an 
ensemble  averaging  technique  [12].  Beard  et  al.  carried  out  a 
combined  experimental  and  computational  investigation  of  the  flow 
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fields  in  a  direct  injection  diesel  engine.  Simulations  were  performed 
with  a  modified  version  of  KIVA  II  code  taking  into  account,  the 
optical  piston  elasticity  due  to  the  high  pressure  [13].  Brunt  and 
Platts  examined  and  quantified  the  main  sources  of  error  in  gross 
heat  release  calculations  using  the  traditional  first  law  heat  release 
model.  The  most  significant  errors  were  due  to  incorrectly  assigned 
ratio  of  specific  heats  and  charge  to  wall  heat  transfer  rates.  An 
alternate  heat  release  model,  the  polytropic  index  first  law  model 
was  developed  and  compared  with  experimental  results  obtained 
from  a  direct  injection  diesel  engine  [14].  Okazaki  et  al.  studied  the 
improvement  of  CFD  approach  for  practical  designing  of  the  direct 
injection  diesel  intake  and  combustion  system.  A  tetrahedron  cells 
model  was  effective  to  get  a  relative  evaluation  of  steady-state  flow. 
The  flow  behavior  of  the  actual  engine  operating  condition  could  be 
predicted  by  the  transient  analysis  that  could  be  used  for 
quantitative  comparison  of  intake  port  performance  [15], 

Bianchi  et  al.  investigated  the  influence  of  different  initial 
condition  procedures  on  combustion  and  emissions  predictions  in 
small-bore  high-speed  direct  injection  diesel  engine.  The  analysis 
was  performed  by  using  the  STAR-CD  code  for  the  intake  stroke 
calculation  and  KIVA-II  for  the  compression  stroke  and  combustion 
simulation.  Analysis  of  the  evolution  of  fluid  dynamic  parameters 
has  revealed  that  a  detailed  definition  of  the  initial  conditions  is 
required  to  properly  predict  the  mean  and  turbulent  flow  field  at 
the  time  of  injection  near  TDC  [16].  Lin  carried  out  numerical 
analysis  of  air  motion  in  three  different  shapes  of  combustion 
chambers  (central  projection  type,  shallow  W  type  and  pataloid 
type)  in  order  to  reduce  the  NO*  and  smoke  emissions.  The 
experimental  study  on  a  single  cylinder  diesel  engine  with  central 
projection  chamber  showed  that  this  system  exhibited  the  rapid 
combustion  process  and  short  main  combustion  period.  Retarda¬ 
tion  of  injection  could  realize  a  large  reduction  in  NO*  and  slight 
rise  in  fuel  consumption  and  smoke  emission  [17]. 

Reitz  et  al.  examined  in  detail  the  velocity  field  in  a  direct 
injection  diesel  engine  employing  both  in-cylinder  velocity 
measurements  and  results  of  multi-dimensional  calculations  using 
KIVA  3V  code  [18].  In-cylinder  fluid  velocity  was  measured  by 
Miles  et  al.  in  an  optically  accessible,  fired  HSDI  engine  at  idle.  The 
velocity  field  was  also  calculated,  including  the  full  induction 
stroke,  using  multi-dimensional  fluid  dynamics  and  combustion 
simulation  models.  A  detailed  comparison  between  the  measured 
and  calculated  velocities  was  performed  to  validate  the  computed 
results  and  to  gain  a  physical  understanding  of  the  flow  evolution. 
Motored  measurements  were  also  presented,  to  clarify  the  effects 
of  the  fuel  injection  process  and  combustion  on  the  velocity  field 
evolution.  Modification  of  the  mean  flow  by  fuel  injection  and 
combustion  was  also  performed.  Substantial  changes  in  the 
tangential  velocity  near  the  spray  path  were  observed  during 
the  fuel  injection  event  and,  after  combustion,  a  radial  profile  of  the 
mean  tangential  velocity  that  favours  turbulence  production  was 
developed  within  the  bowl  [19].  Hountalas  and  Pariotis  developed 
a  simplified  model  for  the  spatial  distribution  of  temperature  in  a 
motored  direct  injection  diesel  engine.  For  the  fluid  flow,  a 
phenomenological  model,  which  was  coupled  to  a  CFD  method  to 
solve  the  energy  conservation  equation  and  therefore  the 
temperature  field,  was  used.  The  proposed  method  has  the 
advantage  of  simplicity  and  low  computational  time.  Experiments 
were  conducted  on  a  direct  injection  diesel  engine  with  bowl  in 
piston  at  various  speeds  and  the  experimental  compression  curve 
was  compared  with  the  theoretical  one.  A  very  good  agreement 
between  the  predicted  and  the  experimental  cylinder  pressure  was 
observed  [20]. 

Dong  et  al.  established  advanced  user  subroutines  and  overall 
simulation  strategies  to  model  a  engine  in-cylinder  turbulent  flow, 
temperature,  pressure,  and  exhaust  gas  recirculation  (EGR) 


concentration  fields  and  to  simulate  EGR  stratification  process 
in  a  typical  pent-roof  gasoline  engine  cylinder  during  intake  and 
compression  strokes.  A  series  of  CFD  studies— combustion  dome, 
ports,  valve  and  valve  mask  configurations,  EGR  tolerances,  swirl 
ratio,  boundary  and  initial  conditions,  etc.,  had  been  conducted  to 
relate  the  “corndog”  flow  structure  and  physical  parameters— 
centrifugal  force,  heat  transfer,  swirl  ratio,  tumble  flow,  turbulent 
and  burnt  gas  residual.  Another  method  to  stratify  EGR  in  cylinder, 
named  “conjugate  vortex  EGR  stratification”  system,  was  dis¬ 
covered  by  the  CFD  analysis.  A  clear  stratification  interface 
between  the  two  streams  of  flow  was  found  and  located  in  the 
cylinder  center  [21].  Auriemma  et  al.  illustrated  an  experimental 
and  numerical  investigation  of  the  flow  generated  by  an  intake  port 
model  for  a  heavy-duty  direct  injection  (HDDI)  diesel  engine.  Tests 
were  carried  out  on  a  steady-state  airflow  test  rig  to  evaluate  the 
global  fluid-dynamic  efficiency  of  the  intake  system,  made  by  a 
swirled  and  a  directed  port,  in  terms  of  mass  flow  rate,  flow 
coefficients  and  swirl  number.  In  addition,  the  LDA  technique  was 
applied  to  obtain  the  local  distribution  of  the  air  velocity  within  a 
test  cylinder.  The  computation  was  carried  out  by  the  fluid- 
dynamic  code  STAR-CD  that  solves  the  ensemble  averaged 
conservation  equations  for  mass,  momentum  and  energy  in 
steady-state  conditions  with  the  turbulence  model.  The  grid, 
reproducing  the  geometry  of  the  intake  port  and  the  real  fluid 
system,  was  made  using  CAD  data  and  the  boundary  conditions 
were  the  same  as  the  experimental  ones.  Results  of  the  computed 
mass  flow  rate,  flow  coefficients,  and  velocity  profiles  was 
compared  to  the  experimental  ones  [22]. 

The  means  of  generating  the  initial  values  and  boundary  data  was 
presented  and  the  effect  of  different  methods  was  discussed  by 
Antila  et  al.  This  study  deals  with  three  different  compression 
ignition  engines  with  cylinder  diameters  of  111,  200  and  460  mm. 
The  initial  cylinder  charge  had  been  carefully  analyzed  through  gas 
exchange  pressure  recordings  and  corresponding  one-dimensional 
simulation.  The  swirl  generated  by  intake  ports  in  a  high-speed 
engine  was  simulated  and  measured.  The  combustion  simulation 
using  a  whole  cylinder  model  was  compared  with  a  sector  model 
simulation  result.  The  difference  between  a  sector  mesh  simulation 
and  a  whole  cylinder  model  simulation  was  found  to  be  considerable 
in  case  of  a  high-speed  diesel  engine.  In  case  of  a  medium-speed 
engine  the  gains  of  using  the  whole  cylinder  model  were  more 
questionable.  The  simulation  of  fuel  injection  rate  and  injection 
velocity  at  the  nozzle  exit  was  presented  in  this  study,  as  well.  The 
validation  of  this  simulation  was  made  through  injection  pressure 
and  injector  needle  lift  measurements.  The  fuel  injection  boundary 
condition,  the  predicted  injection  velocity,  was  found  to  have  a 
noteworthy  effect  on  the  simulated  heat  release  [23].  Time-resolved 
digital  particle  image  velocimetry  (TRDPIV)  data  was  presented  for 
the  in-cylinder  flow  field  of  a  motored  four-stroke,  multi-valve 
direct  injection  spark  ignition  (DISI)  optical  internal  combustion  (IC) 
engine  by  Jarvis  et  al.  The  authors  highlighted  the  application  of  a 
new  TRDPIV  system  to  provide  both  spatial  and  temporal  in-cylinder 
flow  field  development  over  multiple  engine  cycles  for  improved 
understanding  of  cyclic  variation  in  this  paper  [24]. 

Lee  and  Zhao  modeled  the  blow-by  phenomenon  in  a  Small- 
Bore  High-Speed  Direct-Injection  Optically  Accessible  Diesel 
Engine.  They  observed  that  the  difference  in  peak  pressure  was 
up  to  25%  near  TDC.  To  account  for  the  pressure  differences,  a  0-D 
crevice  flow  model  with  a  dynamic  ring  pack  model  was 
incorporated  into  the  KIVA  code  to  determine  the  amount  of 
blow-by.  The  crevice  flow  model  takes  into  consideration  the  flow 
through  the  circumferential  gap,  ring  gap,  and  the  ring  side 
clearance.  Chemical  reactions  of  combustion  were  calculated  by 
using  the  shell  ignition  model  and  KIVA’s  combustion  subroutine. 
Pressure  traces  for  each  case  showed  dramatic  improvement  with 
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the  use  of  the  crevice  flow  model  and  reported  values  around  4%  of 
blow-by  has  seen.  One  of  the  main  causes  for  such  high  blow-by 
values  found  was  attributed  to  the  lack  of  an  oil  film  layer  on  the 
surface  of  the  cylinder  wall  [25]. 

3.  Turbulence  modeling 

Most  of  computational  works  on  diesel  engines  employ  two- 
equation  turbulence  models,  predominantly  the  standard  k-e 
model  for  modeling  turbulence  phenomena.  A  brief  review  has 
been  made  about  these  models. 

A  flow  model  was  presented  by  Borgnakke  et  al.,  which  predicts 
the  swirl  and  turbulent  velocities  in  an  open  chamber,  cup-in- 
piston  IC  engine.  The  swirl  model  was  based  on  an  integral 
formulation  of  the  angular  momentum  equation  solved  with  as 
assumed  tangential  velocity  profile  form,  V\gu(r).  This  enables  the 
swirl  model  to  predict  a  non-solid  body  rotation,  which  is  a 
function  of  the  inlet  flow,  wall  shear  and  squish  motion  during  the 
engine  cycle  [26]. 

A  computer  simulation  for  the  performance  of  a  four-stroke 
spark-ignition  engine  was  used  to  assess  the  effects  of  in-cylinder 
flow  processes  on  engine  efficiency.  The  engine  simulation  model 
is  a  thermodynamic  model  coupled  to  submodels  for  the  various 
physical  processes  of  in-cylinder  swirl,  squish  and  turbulent 
velocities,  heat  transfer  and  flame  propagation.  The  swirl  and 
turbulence  models,  account  for  the  effects  of  changes  in  geometry 
of  the  intake  system  and  the  chamber  design  on  in-cylinder  flow 
processes.  The  combustion  model  was  an  entrainment  burn-up 
model  applicable  to  the  mixing  controlled  region  of  turbulent 
flame  propagation.  The  flame  was  assumed  to  propagate 
spherically  from  one  or  two  spark  plug  locations.  A  heat  transfer 
model  that  was  dependent  upon  the  turbulence  level  was  used  to 
compute  the  heat  loss  from  the  unburned  and  burned  gases.  These 
submodels  are  calibrated  from  experimental  data  by  Davis  and 
Borgnakke  [27].  Musculus  and  Rutland  developed  a  turbulent 
combustion  model  based  on  the  coherent  flame  model  and 
applied  it  to  diesel  engines.  The  combustion  was  modeled  in  three 
distinct  but  overlapping  phases:  low  temperature  ignition 
kinetics  using  the  shell  ignition  model,  high  temperature 
premixed  burn  using  a  single  step  Arrhenius  equation  and  the 
flamelet  based  diffusion  burn.  The  characteristics  of  the  modeling 
approach  were  addressed  by  examining  the  spatial  resolution  of 
the  model,  results  inside  the  engine  cylinder.  The  location  and 
magnitude  of  model  heat  releases,  flame  areas,  and  spray 
distribution  were  examined  [28]. 

McLandress  et  ai.  simulated  a  heavy-duty  caterpillar  3406 
diesel  engine  using  K1VA-3  and  compared  two  turbulence  models 
(k-s  and  RNG  k-s).  The  results  for  the  two  models  did  not  differ 
significantly.  When  the  RNG  k-e  model  was  used,  the  vortex  in  the 
upper  far  corner  away  from  the  inlet  valve  had  a  smaller  radius 
than  when  the  standard  k-e  model  was  used.  The  large-scale  fluid 
motion,  which  developed  during  the  intake  stroke,  was  the  most 
significant  component  of  mixing  during  the  entire  process.  Vortex 
tracking  elucidated  the  survival  of  several  coherent  large-scale 
structures,  which  caused  enhanced  mixing  during  the  latter  half  of 
the  compression  stroke.  These  structure  also  lead  to  a  very 
disordered  and  inhomogeneous  in-cylinder  flow  field  just  prior  to 
the  time  fuel  would  normally  be  injected  into  the  engine.  Two 
mixing  parameters  were  very  helpful  in  characterizing  the  small 
scale  motion  of  turbulence.  Turbulent  kinematic  viscosity  was 
significant  because  it  was  the  primary  indicator  of  spatial  locations 
which  are  favourable  for  the  survival  of  large-scale  coherent 
structure.  The  production  of  turbulent  kinetic  energy  was  also 
helpful  by  capturing  the  locations  where  turbulence  was  being 
created  such  that  small-scale  turbulence  could  be  enhanced  [29]. 


Jakirlic  et  al.  investigated  the  potential  of  the  single-point  closures 
for  predicting  the  flow  and  turbulence  in  piston-cylinder  config¬ 
urations.  The  study  of  Reynolds  Averaged  Navier  Stokes  modeling 
is  simplified  in  simplified  internal  combustion  engine  configura¬ 
tions  [30].  Bianchi  et  al.  proposed  improvements  in  the  RNG  k-e 
model  to  consider  the  effect  of  non-equilibrium  turbulence,  while 
modeling  spray  combustion  in  diesel  engines.  The  simulations 
showed  that  the  two-equation  linear  eddy  viscosity  models  failed 
in  predicting  the  influence  of  turbulence  on  the  fuel  conservation 
rate.  Analysis  of  the  variation  of  the  molecular  viscosity  and 
turbulence  eddy  turnover  time  to  mean  flow  time  scale  ratio 
indicated  that  engine  flows  during  compression,  spray  injection 
and  combustion  undergo  non-equilibrium  conditions.  Three 
approaches  had  been  applied  in  order  to  account  for  non¬ 
equilibrium  turbulence.  The  first  approach  assumes  that  the 
turbulent  time  scale  r  under  non-equilibrium  conditions  is 
proportional  to  the  turbulent  Reynolds  number,  to  account  for 
both  energy  dissipation  from  larger  to  smaller  scales  and  energy 
dissipation  into  heat  at  small  scales.  The  second  approach  assumes 
that  the  turbulent  time  scale  under  non-equilibrium  conditions  is 
in  addition  proportional  to  the  ratio  between  mass  averaged 
turbulent  time  scale  evaluated  under  non-equilibrium  and  mass 
averaged  eddy  turnover  time  evaluated  under  equilibrium.  The 
third  approach  assumes  that  the  turbulent  time  scale  under  non¬ 
equilibrium  conditions  is  proportional  to  the  square  root  of  the 
turbulent  Reynolds  number,  to  project  that  there  occurs  a  delay  in 
response  of  the  turbulent  cascade  to  the  changes  in  the  mean  flow. 
The  first  and  second  approaches  improved  combustion  predictions 
over  a  wide  range  of  operating  conditions  whereas  the  third  was 
not  able  to  show  much  improvement  [31].  A  modified  version  of 
the  laminar  and  turbulent  characteristic  time  combustion  model 
and  the  Hiroyasu-Magnussen  soot  model  was  implemented  in  the 
flow  solver  Star-CD  by  Kaario  et  al.  Combustion  simulations  of 
three  direct  injection  diesel  engines,  utilizing  the  standard  k-\ge 
turbulence  model  and  a  modified  version  of  the  RNG  k-\ge 
turbulence  model,  was  performed  and  evaluated  with  respect  to 
combustion  performance  and  emissions.  Adjustments  of  the 
turbulent  characteristic  combustion  time  coefficient,  which  were 
necessary  to  match  the  experimental  cylinder  peak  pressures  of 
the  different  engines,  had  been  justified  in  terms  of  non¬ 
equilibrium  turbulence  considerations.  The  results  confirm  the 
existence  of  a  correlation  between  the  integral  length  scale  and  the 
turbulent  time  scale.  This  correlation  can  be  used  to  predict  the 
combustion  time  scale  in  different  engines.  It  was  found  that, 
although  the  standard  k-\ge  turbulence  model  produced  adequate 
results  for  the  larger  engines,  the  RNG  turbulence  model  gave 
better  agreement  with  the  experimental  data  for  all  engines  [32]. 

Bianchi  et  al.  focused  on  the  influence  of  the  eddy  viscosity 
turbulence  models  (EVM)  in  CFD  three-dimensional  simulations  of 
steady  turbulent  engine  intake  flows  in  order  to  assess  their 
reliability  in  predicting  the  discharge  coefficient.  The  results  also 
compared  with  experimental  measurements  at  the  steady  flow 
bench.  High  Reynolds  linear  and  non-linear  and  RNG  k-\ge  models 
had  been  used  for  simulation,  revealing  the  strong  influence  of 
both  the  constitutive  relation  and  the  \ge-  equation  formulation  on 
the  obtained  results,  while  limits  in  the  applicability  of  more 
sophisticated  near-wall  approaches  are  briefly  discussed  in  the 
paper  [33].  Miles  et  al.  performed  combined  experimental  and 
computational  approach  was  taken  to  identify  the  source  of  this 
turbulence.  The  results  indicated  that  the  dominant  source  of  the 
increased  turbulence  was  associated  with  the  formation  of  an 
unstable  distribution  of  mean  angular  momentum,  characterized 
by  a  negative  radial  gradient  [34].  Bianchi  and  Fontanesi 
demonstrate  that  RANS  simulation  can  provide  reliable  design- 
oriented  results  by  using  low-Reynolds  cubic  k-\ge  turbulence 
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models.  Different  engine  intake  port  assemblies  and  pressure 
drops  had  been  simulated  by  using  the  CFD  STAR-CD  code. 
Computations  was  performed  by  comparing  two  turbulence 
models:  low-Reynolds  cubic  k-\ge  and  high-Reynolds  cubic  k- 
\ge.  The  analysis  leads  to  remarkable  assessments  in  the  definition 
of  a  correct  and  reliable  methodology  for  the  evaluation  of  engine 
port  breathing  capabilities.  Comparison  between  numerical  results 
and  experiments  showed  that  the  low-Reynolds  cubic  k-\ge  model 
was  unavoidable  to  correctly  capture  the  influence  of  port  feature 
variations  on  engine  permeability  [35]. 

The  one-equation  subgrid  scale  model  for  the  large  Eddy 
simulation  (LES)  turbulence  model  had  been  compared  to  the 
popular  k-\ge  RNG  turbulence  model  in  very  different-sized  direct 
injection  diesel  engines  by  Kaario  et  al.  [36].  Murad  et  al.  compared 
the  suitability  of  various  turbulence  models  in  simulating  the  flow. 
Commercial  softwares  (FLUENT  and  SWIFT)  were  used  to  compare 
the  performance  of  various  turbulence  models.  In  FLUENT,  a 
simplified  vehicle  model  with  slanted  A-pillar  geometry  was 
generated  using  GAMBIT  and  in  SWIFT,  the  simplified  vehicle 
model  was  generated  using  Fame  Hybrid.  CFD  simulations  were 
carried  out  using  FLUENT  under  steady-state  conditions  using 
various  turbulence  models  (k-,  k-Realize,  k-RNG,  k-  and  Spalart 
Allamaras).  In  SWIFT,  k-,  A-RSM  and  HTM2  turbulence  models 
were  used  for  the  steady-state  simulations.  Investigations  were 
carried  out  at  velocities  of  60,  100  and  140  km/h  and  at  0°  yaw 
angle.  Results  of  Cp  values  were  compared  against  available 
experimental  data.  Qualitative  flow  visualization  analyses  based 
on  the  CFD  results  were  also  carried  out  to  better  understand  the 
airflow  behavior  behind  the  A-pillar  region.  Results  obtained  using 
CFD  model  provided  reasonable  agreement  against  available 
experimental  data  [37].  Bai  studied  the  turbulence  field  in  an 
optical  diesel  engine  operated  under  motored  conditions  using 
both  LES  and  particle  image  velocimetry  (PIV).  The  study  was 
performed  in  a  laboratory  optical  diesel  engine  based  on  a  recent 
production  engine  from  VOLVO  Car.  PIV  was  used  to  study  the  flow 
field  in  the  cylinder,  particularly  inside  the  piston  bowl  that  is  also 
optical  accessible.  LES  was  used  to  investigate  in  detail  the 
structure  of  the  turbulence,  the  vortex  cores,  and  the  temperature 
field  in  the  entire  engine,  all  within  a  single  engine  cycle.  The  LES 
results  are  compared  with  the  PIV  measurements  in  a 
40  mm  x  28  mm  domain  ranging  from  the  nozzle  tip  to  the 
cylinder  wall.  In  the  intake  phase  the  large-scale  swirling  and 
tumbling  flow  streams  are  shown  to  be  responsible  for  the 
generation  of  large-scale  vortex  pipes,  which  break  down  to  small- 
scale,  turbulent  eddies.  In  the  later  phase  of  compression, 
turbulence  was  mainly  produced  in  the  engine  bowl.  Turbulence 
and  the  large-scale  coherent  vortex  shedding  due  to  the  Kelvin- 
Helmholtz  instability  are  responsible  for  the  enhanced  heat  transfer 
between  the  bulk  flow  and  the  walls.  A  temperature  in  homogeneity 
of  about  50-60  K  can  be  generated  in  the  cylinder  [38]. 

PIV  and  paddle  wheel  measurements  have  been  conducted  in  a 
static  heavy-duty  diesel  engine  rig  to  characterize  the  flow  features 
with  different  valve  lifts  and  pressure  differences  by  Kaario  et  al. 
These  measurements  were  compared  with  CFD  predictions  of  the 
same  engine.  The  simulations  were  done  with  the  standard  k-\ge 
turbulence  model  and  with  the  RNG  turbulence  model  using  the 
Star-CD  flow  solver  [39].  Lendormy  et  al.  studied  initial  conditions 
for  turbulence  prior  to  combustion  with  the  help  of  a  four-valve, 
large  bore  diesel  engine  CFD  model.  A  one-dimensional  model  was 
first  used  to  provide  boundary  conditions  as  well  as  the  initial  flow 
conditions  at  the  beginning  of  the  simulation.  Steady-state  and 
transient  boundary  conditions  were  studied.  In  addition,  the 
standard  k-\ge  and  RNG/k-\ge  turbulence  models  were  compared. 
From  the  averaged  values  of  turbulence  kinetic  energy  and  its 
dissipation  rate  over  the  cylinder  volume,  a  re-tuned  correlation 


for  defining  the  initial  turbulent  conditions  at  bottom  dead  center 
(BDC)  prior  to  the  compression  stroke  was  proposed  [40]. 

4.  Spray  modeling 

The  simulation  of  evaporating  droplet  dynamics  in  diesel 
engine  simulation  are  based  mostly  on  models  proposed  for  single 
droplet  evaporation  phenomena.  A  brief  review  of  spray  model  has 
been  made. 

Haselman  and  Westbrook  developed  a  theoretical  model  for 
calculations  of  the  evolution  of  fuel  injected  into  internal 
combustion  engine  chambers.  Fuel  injection  in  the  form  of  a 
gaseous  jet  and  in  the  form  of  a  liquid  droplet  spray  is  considered. 
The  model  used  the  basic  equations  of  conservation  of  mass, 
momentum,  and  energy  in  both  the  gaseous  and  liquid  phases. 
Applications  in  two-dimensional  symmetry  of  the  gaseous  jet  form 
of  the  numerical  model  are  described  in  stagnant  and  swirling 
airflows.  The  liquid  droplet  spray  model,  including  coupling 
between  the  droplets  and  the  gas  phase  medium,  are  described. 
Applications  of  the  two-phase  model  are  described  for  the  case  of 
axially  symmetric  injection.  Finally,  the  liquid  spray  model  and  the 
gas  jet  model  are  applied  to  the  same  conditions,  leading  to  a 
general  assessment  of  the  ranges  of  validity  of  the  gas  jet  model 
[41]- 

An  analytical  spray  model  was  described  which  can  be  used  to 
calculate  the  penetration  and  trajectory  of  a  spray  in  an  engine 
combustion  chamber  with  air  swirl  by  Sinnamon  et  al.  The  model 
consists  of  integral  continuity  and  momentum  equations  written 
for  a  steady-state  gas  jet.  The  model  contains  adjustable 
entrainment  and  drag  parameters  evaluated  from  experimental 
data.  A  special  single-cylinder,  see-through  engine  and  a  schlieren 
optical  system  were  used  to  study  transient  liquid  fuel  sprays 
under  varied  conditions.  These  experimental  observations  were 
used  to  determine  appropriate  values  for  the  adjustable  para¬ 
meters  in  the  spray  model.  Comparisons  between  model  calcula¬ 
tions  and  the  experimentally  observed  sprays  were  presented  for  a 
wide  range  of  conditions  [42].  Liu  et  al.  modeled  the  effects  of  drop 
drag  and  breakup  on  fuel  sprays.  Spray  atomization  was  modeled 
using  a  wave  breakup  model  that  was  based  on  results  from  jet 
stability  analysis,  and  also  using  KIVA’ s  Taylor  analogy  breakup 
(TAB)  model.  The  wave  model  was  found  to  give  good  results  for 
both  drop  trajectories  and  breakup  drop  sizes.  The  best  results 
were  obtained  with  a  smaller  value  of  the  breakup  time  model 
constant  than  previously  used  in  spray  computations.  Drop 
trajectory  and  size  measurement,  together  with  high  magnifica¬ 
tion  photographs,  indicated  that  drop  distortion  should  be 
accounted  for  in  sprays.  Accordingly,  a  modified  drop  drag  model 
was  proposed  in  which  the  drop  drag  coefficient  changed 
dynamically  with  the  flow  conditions  during  the  drop  lifetime. 
In  the  model,  the  value  of  the  drag  coefficient  varied  between  the 
limits  of  a  rigid  sphere  (no  distortion)  and  a  disk  (maximum 
distortion).  The  drop  distortion  was  computed  using  the  Taylor 
analogy  between  a  drop  and  a  spring-mass  system,  and  the 
breakup  process  was  described  using  the  wave  model.  The  TAB 
breakup  model  was  also  found  to  give  good  predictions  of  drop 
trajectories,  but  the  model  under  predicted  measured  the  breakup 
drop  size  considerably.  The  TAB  results  were  thus  relatively 
insensitive  to  the  drop  drag  model,  since  small  drops  have  low 
inertia,  and  they  quickly  accelerated  up  to  the  gas  velocity.  The 
wave  breakup  and  dynamic  drag  models  were  also  applied  to 
diesel  sprays.  The  results  confirmed  previous  studies  that  show 
that  spray-tip  penetration  is  relatively  insensitive  to  drop  breakup 
and  drag  models.  However,  the  distribution  of  drop  sizes  within 
the  sprays  was  found  to  be  influenced  by  the  model  details.  This 
was  due  to  the  fact  that  drop  drag  changes  the  drop-gas  relative 
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velocity,  and  this  changes  the  spray  drop  size,  since  the  drop 
breakup  and  coalescene  processes  depend  on  the  velocity. 
However,  these  changes  occurred  in  such  a  way  that  the  net 
effect  on  the  spray  penetration  is  small  [43]. 

Senda  et  al.  proposed  new  submodels  for  the  dispersion  process 
of  the  spray  impinging  on  a  flat  wall  to  reproduce  the  spray/wall 
interaction  process  inside  the  diesel  combustion  chamber.  These 
models  treated  the  non-vaporizing  spray  impinging  on  the  wall  at 
room  temperature  and  the  vaporizing  spray,  which  impinges  on 
the  wall  at  a  high-temperature  above  the  spray  saturation 
temperature.  These  models  were  incorporated  into  the  KIVA  code 
to  calculate  the  dispersion  process  [44].  Bai  and  Gosman  developed 
and  assessed  a  spray  impingement  model  for  diesel  spray 
simulation.  Validation  of  model  had  been  performed  against 
several  sources  of  experimental  data.  The  conditions  investigated 
involved  ambient  and  elevated  pressures  and  temperatures,  and 
normal  and  oblique  impingement  cases,  as  well  as  an  engine.  In  all 
the  cases,  the  predictions  compared  reasonably  well  with  the 
experiments,  and  in  some  cases  very  good  agreement  was 
achieved.  They  concluded  that  the  model  presented  was  useful 
in  predicting  spray  impingement  events  [45].  Zhang  et  al. 
introduced  a  very  simple  method  to  measure  the  liquid  phase  of 
spray  in  an  optically  accessible  direct  injection  diesel  engine.  By 
using  this  method,  the  effects  such  as  injection  pressure,  nozzle 
specification,  intake  air  boost  and  temperature  on  liquid  phase 
penetration  before  ignition  was  investigated  [46].  Pelloni  and 
Bianchi  proposed  a  hybrid  model  for  simulation  of  high  pressure 
driven  diesel  fuel  spray.  The  extensive  comparison  between 
numerical  results  and  experimental  has  shown  the  reasonably 
good  accuracy  of  the  hybrid  model  in  predicting  the  breakup 
process  of  high  density,  high  pressure  driven  sprays.  Poor 
prediction  of  spray  penetration  had  been  detected  by  using  a  grid 
resolution  of  practical  internal  combustion  engine  computations. 
An  anisotropic  model  of  the  droplet  turbulent  dispersion  in 
conjunction  with  the  k-s  model  had  been  presented.  Significant 
changes  in  the  spray  structure  and  droplet  size  were  detected  by 
comparing  the  isotropic  and  the  anisotropic  approaches,  especially 
near  the  nozzle  exit  [47]. 

Bensler  et  al.  performed  simulations  modeling  the  spray 
injection  of  a  5-hole  diesel  injector  in  a  pressure  chamber.  A 
graphical  methodology  was  utilized  to  match  the  spray  resulting 
from  the  widely  used  discrete  droplet  spray  model  to  pressure 
chamber  spray  images.  Using  engine  measurements  to  initialize 
the  combustion  chamber  conditions,  the  compression  stroke,  the 
spray  injection  and  the  combustion  simulation  was  performed.  The 
novel  RTZF  two-zone  flamelet  combustion  model  was  used  for  the 
combustion  simulation  and  was  tested  for  partial  load  operating 
conditions  [48]. 

Montajir  et  al.  investigated  the  fuel  spray  development  under 
high  injection  pressure  and  high  gas  charging  pressure  photo¬ 
graphically  in  a  small  direct  injection  diesel  engine  with  a  common 
rail  injection  system.  Improvement  in  the  spray  distribution  inside 
the  cavity  was  attempted  by  changing  the  lip  shape,  side  wall 
shape,  and  injection  angle  [49].  Bae  et  al.  investigated  the  spray 
characteristics  from  diesel  injectors  with  various  nozzle  geome¬ 
tries  by  microscopic  and  macroscopic  visualization.  Spray  angle 
measured  from  the  macroscopic  images  were  hardly  affected  by 
the  injection  pressure  at  the  atmospheric  conditions.  The  spray 
angle  measured  from  the  microscopic  images  near  the  nozzle  was 
also  independent  of  injection  pressure  [50]. 

Moriyoshi  et  al.  experimentally  evaluated  spray  model  for  a 
swirl-type  injector.  Volume  of  fluid  model  (VOF)  was  used  to 
simulate  the  two-phase  flow  inside  the  injector  and  also  the  liquid 
film  formation  process  outside  the  nozzle,  while  discrete  droplet 
model  was  used  to  simulate  a  free  fuel  spray  in  a  constant  volume 


chamber  using  initial  conditions  deduced  by  empirical  equations 
or  calculated  results  of  VOF  model  [51].  Chryssakis  and  Assanis 
developed  a  comprehensive  model  for  sprays  emerging  from  high- 
pressure  swirl  injectors  in  D1SI  engines  accounting  for  both 
primary  and  secondary  atomization.  The  model  considers  the 
transient  behavior  of  the  pre-spray  and  the  steady-state  behavior 
of  the  main  spray.  The  pre-spray  modeling  was  based  on  an 
empirical  solid  cone  approach  with  varying  cone  angle.  The  main 
spray  modeling  was  based  on  the  liquid  instability  sheet 
atomization  (LISA)  approach,  which  are  extended  here  to  include 
the  effects  of  swirl  [52].  An  effort  had  been  made  to  understand  the 
physical  phenomenon  behind  combustion  in  direct  injection  diesel 
engine  and  to  evaluate  its  performance  through  a  mathematical 
model  by  Udayakumar  and  Anand.  The  model  developed  consists 
of  two  phases— the  spray  model  and  the  combustion  model.  In  the 
spray  model  a  2D  axis-symmetric  grid  was  considered.  The  spray 
dispersion,  penetration,  and  the  cells  that  are  ready  for  ignition  at 
the  end  of  ignition  delay  are  found  from  the  spray  model.  A 
numerical  technique  is  used  to  solve  the  conservation  equations  of 
mass,  momentum  and  energy  in  the  combustion  model  to  evaluate 
the  properties  at  each  crank  angle.  The  model  was  also  used  to 
determine  mass  fractions  burnt,  pressure,  temperature,  heat 
release  at  each  crank  angle  and  performance  parameters.  To 
validate  the  simulation  results,  test  had  been  conducted  on  a 
single-cylinder,  four-stroke  Kirloskar  direct  injection  diesel  engine 
and  a  reasonable  agreement  was  found  between  the  computed  and 
the  experimental  results  [53].  Strauss  and  Zeng  found  that 
different  engines  could  exhibit  significantly  different  behaviors 
with  similar  fuel  sprays.  This  difference  is  attributed  to  the 
difference  in  scavenging  flow  patterns  and  its  effect  on  the 
momentum  balance  between  the  fuel  spray  and  the  airflow.  Two 
different  sprays  were  tested  in  the  fixture  and  in  a  variety  of 
engines  [54]. 

Lucchini  et  al.  developed  a  comprehensive  model  for  sprays 
emerging  from  high-pressure  swirl  injectors  for  gasoline  direct 
injection  (GDI)  engine.  The  primary  and  secondary  atomization 
mechanism  as  well  as  the  evaporation  process  both  in  standard 
and  superheated  conditions  are  taken  into  account.  The  spray 
modeling  after  the  injection  was  based  on  the  LISA  approach, 
modified  to  correctly  predict  the  liquid  sheet  thickness  at  the 
breakup  length.  The  effect  of  different  values  of  the  superheat 
degree  on  evaporation  and  impact  on  the  spray  distribution  and 
fuel-air  mixing  was  also  analyzed  [55].  Yi  and  DeMinco  developed 
an  engine  simulation  tool  in  a  commercial  CFD  code  to  study  the 
spray  and  mixing  process  that  can  be  used  to  access  the 
performance  of  a  GDI  engine  [56].  An  experimental  and  numerical 
study  is  conducted  on  the  spray  and  mixture  properties  of  a  hole- 
type  injector  for  direct  injection  gasoline  engines.  The  laser 
absorption  scattering  (LAS)  technique  was  adopted  to  simulta¬ 
neously  measure  the  spatial  concentration  distributions  and  the 
mass  of  the  liquid  and  vapor  phases  in  the  fuel  spray  injected  into  a 
high-pressure  and  high-temperature  constant  volume  vessel.  The 
experimental  results  are  compared  to  the  numerical  calculation 
results  using  three-dimensional  CFD  and  the  multi-objective 
optimization.  In  the  numerical  simulation,  the  design  variable  of 
the  spray  model  was  optimized  by  choosing  spray-tip  penetration, 
and  mass  of  liquid  and  vapor  phases  as  objective  functions  by  Sato 
et  al.  [57].  An  experimental  study  of  real  multi-hole  diesel  nozzles 
with  different  geometry  in  terms  of  conicity  factor  was  performed 
under  current  direct  injection  diesel  engines  operating  conditions 
by  Payri  et  al.  A  complete  characterization  of  the  internal  nozzle 
flow  and  the  sprays  injected  under  non-evaporative  conditions  had 
been  performed  [58].  Xue  et  al.  developed  a  methodology  that  uses 
a  locally  refined  mesh  in  the  spray  region.  This  methodology  used 
an  h-refinement  adaptive  method.  In  addition,  they  demonstrated 
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the  efficiency  and  accuracy  of  the  adaptive  mesh  refinement 
scheme  [59].  Abani  et  al.  proposed  an  improved  spray  model  for 
reducing  numerical  parameter  dependencies  in  diesel  engine  CFD 
simulations.  The  improved  spray  model  was  implemented  in  the 
engine  simulation  code  KIVA-3V  and  was  tested  under  various 
conditions,  including  constant  volume  chambers  and  various 
engine  geometries  with  vaporizing  and  combusting  sprays  with 
detailed  chemistry  [60].  Micklow  and  Abdel-Salam  presented  a 
computationally  efficient  numerical  model  that  accurately  pre¬ 
dicts  complex  spray  distribution  and  spray  penetration  for  a  direct 
injection  compression  ignition  engine  configuration.  A  modified 
version  of  a  multi-dimensional  computer  code  KIVA-3V  is  used  for 
the  computations,  with  improved  submodels  for  varying  mean 
droplet  diameter,  varying  injection  velocity  and  drop  distortion 
and  drag  [61].  Kim  et  al.  investigated  the  air/fuel  mixture 
formation  and  combustion  characteristics  in  a  spray-guided 
GDI  engine  using  a  commercial  code,  STAR-CD.  The  spray  injection 
has  modeled  arranging  multiple  points  using  random  function 
along  the  ring-shaped  nozzle  exit.  The  present  approach  reason¬ 
ably  predicted  the  spray  behavior,  the  mixture  distribution  near 
the  spark  plug,  and  flame  propagation  inside  the  combustion 
chamber  [62]. 

5.  Conclusions 

Improved  computational  mesh  generation  techniques  and 
efficiency  of  solver  have  widely  influenced  the  application  of 
CFD  methods  to  reciprocating  engine  models,  which  require 
movable  domain  boundaries  and  compressible/expandable 
meshes.  Multi-dimensional  models  are  able  to  predict  in-cylinder 
mean  gas  velocities  with  high  accuracy,  but  RMS  components  are 
predicted  with  lesser  accuracy.  The  in-cylinder  flow  measure¬ 
ments  are  performed  with  non-intrusive  methods  such  as  LDA,  PIV 
and  optical  methods,  thus  increasing  the  reliability  of  experi¬ 
mental  data  obtained  inside  the  engine  cylinder.  Most  computa¬ 
tional  works  on  diesel  engines  employ  two-equation  turbulence 
models,  predominantly  the  standard  k-e  model  for  modeling 
turbulence  phenomena.  The  simulation  of  evaporating  droplet 
dynamics  in  diesel  engine  simulations  is  based  mostly  on  models 
proposed  for  single  droplet  evaporation  phenomena.  Ignition 
models  used  are  based  on  the  shell  ignition  model.  A  large  number 
of  combustion  models  are  constituted  for  diesel  combustion 
simulation,  but  do  not  find  further  validation.  The  laminar  and 
turbulent  characteristics  time  scale  model  finds  higher  usage  when 
compared  to  the  other  models  for  simulating  diesel  combustion. 

Most  combustion  simulations  are  concentrated  in  the  mid  and 
high  load  range.  However,  given  the  importance  of  higher  emission 
weighting  factors  at  idling  conditions  in  various  duty  cycle,  and 
until  recent,  the  availability  of  considerable  amount  of  in-cylinder 
experimental  data  obtained  under  firing  conditions  at  low  loads, 
work  on  engine  simulations  are  not  available  at  low  speeds. 
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